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SUMMARY 


The propagation characteristics of a crack in a brittle, linear elastic material have been 
investigated. In order to obtain a mathematical description, the material is supposed to be 
homogeneous and isotropic both as regards stress-strain relations and fracturing characteristics. 
It is & priori assumed that a two-dimensional crack can be propagated and that this crack 
propagates in one plane, only. Further, the surface energy of the material is supposed to be 
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zero. Then the crack-tip velocity will be constant (cf. Mott, 1948), and the problem of the 
crack propagation can be investigated by a study of displacements and stresses in a semi- 
infinite solid, where the plane surface is subjected 


1. to a constant pressure, acting on an infinite strip, the width of which is symmetrically 
increasing from zero with a constant velocity, and 
2. to a pressure outside this strip such that the normal displacement of the surface outside 


this strip is zero. 


This mixed boundary value problem has been treated. The stresses in the solid and the 
normal displacement of the surface have been solved exactly. 

The result shows that the displacement of the surface is elliptic, just as in the corresponding 
static case. By studying the stresses in the vicinity of the edge of the strip and using some 
simple fracture criteria it can be shown that the assumption a priori regarding the existence 
of a two-dimensional crack propagating in one plane is justified. Further it can be concluded 
that the crack-tip velocity equals the velocity of surface (Rayleigh) waves. 

In a real material, the existence of a surface energy does not seem to have much influence 
on the crack propagation except for the very start (Mott, 1948).. The inhomogeneity and the 
anisotropy of the fracturing characteristics will cause deviations from a plane propagation of 
the crack. This fact will reduce the mean velocity of the crack in the direction of propaga- 
tion. The effects of non-brittleness or deviations from Hooke’s law are not discussed. 


Introduction 


By using dimensional analysis Mott (1948) was able to show that the maximum 
propagation velocity of a crack is independent of the surface energy. If we con- 
sider a two-dimensional crack of length 2 a, propagating symmetrically with con- 
stant crack-tip velocity 

da _ 

dives’ 
then the surface energy is proportional to a, whereas the release of stress energy 
due to the propagation of the crack and the kinetic energy resulting from this 
propagation are both proportional to a. Thus the energy balance is simply a 
balance between the released stress energy and the kinetic energy when aco. 
Roberts and Wells (1954) assumed that the stress state in the material sur- 
rounding the crack is approximately equal to the corresponding stress state in the 
static case. Using Mott’s criterion for the energy balance they arrived at the 
conclusion that the maximum crack velocity in a brittle material is about 38 
per cent of the sound velocity. The approximation used by Roberts and Wells 
does not, however, seem to be permissible. 

A two-dimensional crack will be considered in the following. The material is 
assumed to be homogeneous and isotropic both as regards fracturing char- 
acteristics and stress-strain relations. It is supposed to obey Hooke’s law and 
to be perfectly brittle. Further the surface energy is assumed to be zero. This 
assumption simplifies the study, because the crack can be supposed to nucleate 
from an infinitesimally small micro-crack and to propagate with maximum velo- 
city from the start. On the other hand, as the maximum velocity obtained is 
independent of the surface energy this assumption does not violate our possibility 
to find the maximum crack velocity in a brittle material. 
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Fig. 1. Coordinate system and position of crack at time ¢. The negative y-axis points towards 
the reader. 


Consider an infinite medium with the characteristics specified above. Cartesian 
coordinates x, y, 2 are introduced (see Fig. 1). The crack is supposed to nucleate 
from an infinitesimally small micro-crack situated along the y-axis from y= — co 
to y=+co and to propagate symmetrically in positive and negative x-direc- 
tions with the constant crack-tip velocity v, when an applied tensile stress 
reaches the nominal value 


02> do- 


Here we assume a state of plane strain, so that nominally 


v : , 
z= y= do (v= Poisson’s ratio) 


oo 
Try = Tyz = Ter = 9. 


We intend to find the crack-tip velocity v,, the shape of the crack, and the 
stresses in the solid. 


Method of treating the problem 


When treating the problem mathematically it is convenient to superpose a 
state of stress 
0z= — % 
y 


ve eng BaP 


% 
Try = Tyz = Tez =9 


defining a state of plane strain. Inside the crack, where the z-component, only, 
is of importance, we may suppose that the superposition is achieved by the 
action of a uniform liquid pressure q. 
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Fig. 2. Load and surface displacement at time ¢ in Problem A. 


After this state of stress has been superposed, the stresses at infinity are all 
zero. A cut is made along the plane z=0 and the semi-infinite solid z>0 is 
considered. Thus the following problem is created (see Fig. 2): 


Problem A.—A pressure qd, is acting on the infinite strip |x|<vot of the sur- 
face z=0 of a semi-infinite solid z>0. The vertical displacement wu, at the sur- 
face z=0 is zero for |a|>v)t. Find the normal stress o, at the surface z=0 
for |z|>vpt and the displacement wu, at the surface z=0 for |a|<vpt. 

After solving this problem, the original problem can be treated after in- 
vestigation of the stresses in the vicinity of the crack-tip, using a suitable frac- 
ture criterion. 

In order to solve Problem A, the following problem is first solved (see Fig. 3): 


Problem B.—A pressure q is acting on the infinite strip |z|<vt of the sur- 
face z=0 of a semi-infinite solid z>0, whereas no pressure acts on the surface 
z=0 for |z|>vt. Find the displacement w, at the surface z=0. ; 

After solving this problem, a superposition of different pressures q with dif- 
ferent propagation velocities v can be made in such a way that the resulting 
pressure is gy for |~|<v jt and the resulting displacement u, at the surface 
z=0 is zero for |x|>vgt. This leads to an integral equation for the normal 
stress og, on the surface z=0 for |a|>v )¢. Solving this integral equation, the 
first part of Problem A is solved. This part of the problem will be called 
Problem C. 

After Problem C is solved, the second part of Problem A is solved, viz. to 
find the displacement wu, at the surface z=0 for |x|<v)t. This part of the 
problem will be called Problem D. 
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Fig. 3. Load and surface displacement at time ¢ in Problem B. 


In fact, instead of working out the displacement uw, in Problem B, the ac- 
celeration @°u,/dt? will be calculated and used for the superposition in Prob- 
lem C, because this simplifies the calculations. It is known that if é?u,/et?=0 
for |x|>vpt, then u,=0 for |x|>vpt because u,=0 for |a|>cat (ca= propaga- 
tion velocity of irrotational waves). 


Symbols 


xX, y, z: Cartesian coordinates. The semi-infinite solid occupies the space z> 0. 

x, ys z: Unit vectors in positive w-, y- and z-directions, respectively. 

“= | Uz, Uy, Uz|: Displacement vector. Conditions of plane strain are prevailing, so that uy=0 
and uw, and wz are functions of x and z, only. 

Oz, Oy, Fz, Try; Tyz, Tzz: Components of stress tensor. As the solid obeys Hooke’s law, Try =Tyz = 0. 

cq: Propagation velocity of irrotational waves. 

c;: Propagation velocity of equivoluminal waves. 

cs: Propagation velocity of surface (Rayleigh) waves. 


Ca 
t: Time 
T=Cqt 
p: Laplace transform variable with respect to T. 


0 for t<0 
U (t)= 4 4/, for t=0 ¢: Unit step function 
1 for t>0 
d : 
6 (t) Fd (t): Dirac’s delta function. 
$= div a 


p= rot i 
o: Density of the solid. 


Other symbols are defined in the text. 
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Problem B 


A pressure g is acting on the infinite strip |x|<vt of the surface z=0 of a 
semi-infinite solid z>0, whereas no pressure acts on the surface z=0 for |a|>ve. 


Find the displacement u, at the surface z=0. 


The equation of motion in an isotropic elastic solid, 


Mathematical statement of the problem 


2- 


OT 


can be split up into the equations 


On the boundary z=0 the following relation holds 


div grad = 


ig 
— rot rot P= 72° 


ad 


ov 


a) 


Cig grad 6 —k* rot 


7 


C5 = — (1-29) 644+ pxé-28 


ou 


az 


where G=stress vector referring to the surface z=0, i.e. 


where 


Pld 2 sOr wel ie 
o= 
0 for |a|> Br, 


Now, because @ is a function of w and z, only: 


b= ¢ (2, 2) 
P=y (a, 2)°G. 


Thus Eqs. (1), (2), (3), (4) turn into 


164 


2- 


7) 


Ag 


—t- grad $1 


_&¢ 


ioe 


Ox 


oy, OY, 
az. =u?) 


> 


(1) 


(2) 


(4) 


ARKIV FOR FYSIK. Bd 18 nr 10 


1 ey 

Mia = 

yp ka 

6 ou ou 
poe eA L + kts | 2 ie 2—*! &. 
od ja-2 Ht * E =| 4 


Introducing the Laplace transforms 


D=p fe" ddr, W=p [e?' wdr, b= p fe? udt 
0 0 0 


eo 


and S,=p f e **g,dt one obtains the equations 
0 


(S,)e-0 = —qe P78 (5) 

Wr ee i 2) Q (6) 
A®D=yp’® (7) 

= aoe (8) 

qe Mele — (y — DAY) +2 4s = ago ea (9) 
rz=0. (10) 


Solution of Eqs. (6)—(8) 


Excluding waves propagating in negative z-direction, the following general 
solutions of Eqs. (7) and (8) are obtained: 


+00 LrLe 

0 = f Ae-VPt+# eos aada (11) 
+0 = ae 

w= [ Be Et sin axda, (12) 


where A and B are functions of p and «, only. 

After insertion of the general solution in Eqs. (9) and (10), and making use 
of the inversion theorem for Fourier integrals, one obtains the following equa- 
tions for determination of A and B: 
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2 2 3 
(Hr e)atea| Brae - q ey 


2k? 2nkocy peta B 
_ Seay st p* 
aVp+ePAth (Sate) B=0. 


Solution of these equations yields 


pee 2 
* Gunie: Sh Bn eam : amie (13) 
a 2 2 2 2 22 2 
Qnk? oc p +a? B (+a) = Vp YE + ot 
Lei otaas 
i nak 
3 Aes Pp 
eR a i apa 3 : (14) 
(F.+2") =o Vp? + & VB +a! 
From Eq. (6) we get for z=0: 
+00 
p?w,= | [—Vp? +0 A—k? «B] cos exda= 
q 2° p Vp? + 
0 (Fat) — a Vp YE + ot 
Confining ourselves to real, positive values of p, we get after putting 
ta|x|= —pd 
the following expression 
—too 
2 =e qd B 4 e —pe 
p" Wz Seite lar p| ans Ps r(B)e dé, (15) 
4 me eS 


@? 

Ble \eo 

ais pe lads Jal RY feo) Sy ee 
(r2-3) +ayi-% Poe 


The integrand of the integral in Eq. (15) has the branch points @= +|a| and 
0=+|«|/k and four simple poles, @=+|2|/B and 6= +6), where 6,>|a|/k. 
Thus the path of integration may be deformed to the path shown in Fig. 4. 


where 13 ( 
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idm 6 


Fig. 4. Path of integration in the §-plane. 


Thus we obtain 


co 92 


2 Ate Jers on B i —pé we 
pw; Sako da] p| 5 Re F(%3)¢ dé 
0 La a) 


be {lb Pel Fe 
—Re p-ni-Fle(4) ev, 


where the integral is to be taken in the sense of the Cauchy principal value, 
and the proper branch of the function F is found by studying Fig. 4. Inverting 
the Laplace-transform, we get for t>0: 


0 uz q B Z| sya 1 | «| 
— = — Re (= | +=" — = Seed RS 
ox 22k* 0 c4|a| a : (S)+5 = p im F (5 0(p ) cot 
Lice 2 
Integration of this expression in order to find uw, is not necessary for solution 
of Problem C. 


The method used for inverting the Laplace transform was introduced by 
Cagniard (1939). 


Problem C 


A pressure g, is acting on the infinite strip |z|<vgt of the surface z=0 of 
a semi-infinite solid z>0. The displacement wu, at the surface z=0 is zero for 
|x|>v)t. Find the normal stress o, at the surface z=0 for |a|> vt. 


Mathematical statement of the problem 


Tf the stress on the surface is 


0 for |2|<fyt 
do,={ —q (f)df for Byr<|x|<Pr 
0 for Bt <|2|. 
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Eq. (16) gives for t>0 


au a (AaB BB | ner()+ 


Ee Qa k* oc%|x t 
Q a | | lag a f= 03 


+3l2l5(p_lel)g Im F (5) - 3 o(6—-l) 0, Im F (=, : 


Superposing such stresses, do,, and the stress 


_ [= % for |x| <fyt 
=2 10 enters r=|at 


on the surface, in such a way that @u,/d1*=0 for |%|> fot, we get the in- 
tegral equation 


iksgstaernth 


fale Ee = p-Im F (5): Lela ( (B) dB + 


2 
+ tho Re F (5) -0 forebat eles 
Putting 2°/t*=r, B’=s and q' (8)=(s), we obtain for pi<r<1 


1 
Im Ff (*) 
r 


2 ete | Pe ) (r= Bo) Wess Bo (r 2) 4s=2 458, 


oC ent Vs 


or 


1 
Im F (*) m= 
%, 


_ bie aries 1 | (s—fo) 9 (s) uP La 2408 
ne r(!) aren | Pomme Ons “sg i fam s)ds=—2", (17) 


Bo B 


Thus a singular integral equation 
lay 
a(r) p (r)+— | EAU (18) 
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has been obtained, where 


1m F(’) 0 for po<r<k? 
a (r) = ————_~_: (r — 3) = a Soiree 
1 4° V(l—r) (r_-k? < 
Re F(*) a dr B2) for Ber<l 
b(r)=r— fo 


1 
) =3 
Kee Po, 1 Vs Lip ey 
IU MLA s 
Bo 


Integral equations of similar type have been studied by Carleman (1922). 
Further studies are reported in the monograph written by Muskhelishvili (1953). 


Solution of Eq. (18) 
The general solution of Eq. (18) can be written (cf. Muskhelishvili (1953)). 


1 


a (r) bat (s — fo) e? b (s) 
=> K ee =  ——_ — K e d 
le a* (r) + B* (r) xt (r — B2) Va? (r) + B? (r) | (s—r) Va? (s) +b? (s) 4 
B- 
7 (19) 


ee 


(r — B2) Va? (r) +B? (nr) 


shoul 


where JL is an arbitrary constant and 


w(r) has been calculated (see Appendix 1) and one obtains 


eee A ea for Bo<r<k? 
oo} Tbs 7-2-4 Vr) (2) 
— 2 , 
tee r (r— ke) for Kk? <r<l1 
r—Bo V(r—2k*)*+ 16k (1 —1) (r—#’) 


The next step is to calculate the integral in Eq. (19). One obtains (see Ap- 
pendix 2): 
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> (5—B3) e? B(8) 
les Va? (s) + 6 (s) 
Bo 
a (M —r) for po<r<Fk? 

r (r —K2) (1—r) Vir — &°) (r — Bo) 
m(M +m 4 ore 61 —n) (rH) 


for k?<r<l, 


where M is a constant that is not necessary to calculate. 
Thus one obtains the following solution of Eq. (18): 
For po<r<k’: 


(M —r) [(r—2k?)? 443 Vie 


‘nines Kis a 
ia r (rk) V(L=1) (r— Be)? 


S _4hVO—n (2?) 
KG 2k)? 4k V(1l—r) (—r) 


r (r —k2) V(l—1) (r — B3)° 


For k?<r<l: 


(Misr) ise ENG ces (r—2 k?)? 


aiaind) i : 
a r (r— ke) V(1—r) (r — B3) r (r — ks) V(1—1) (r — Bo)” 


This can be written in the following way: 
For po<r<k*: 


(r —2?)? 418 V1 —1) (k?— 1) (r —2%)? 42° V(1—r) (k®—1) 
a yr ‘ rt A! 
Ad ae r V(l—r) (r—B3) te (r— k) V(L—1) (r — 3)? a 


For k?<r<l: 


(r —2 k?)? iarvee (r—2 k*)? 


a G09 BF 


if vp V(l—r) (r — B5)° 


(22) 


where Q and Q, are constants. 


It is possible to show (see Problem D and Appendix 3) that the problem is . 
satisfied by the term containing Q, only. Thus the solution is 


(r—2hP—418 V(1—7) (=r) 


Q ——— for Bo<r<k (23) 
rV(1—r) (r— BR) 
gy (r)= ae 
f wz 2 
2 Va—n (r—A) i) a al a (24) 


where Q is to be determined from the equation 
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s)ds 


1 * Va— By 2% Bo. 
~{ Vs Pp ( — 1h 


One obtains (see Appendix 4): 


where 


g (Bo) = (1 —4 2) 3+ 494] K (V1 Ba) — pal 412 (1+ i) 62484] E (VIA) — 
0 


—418 «1~ iy w(Y/ 1 (2) (22)")4 +Sra- 63) B (/1-(&)), (25) 


where K and E are the complete elliptic integrals of the first and the second 
kind, respectively. 


The stress oz (x) at the surface z=0 


The stress o,(x) on the surface z=0 for |x|> fort is given by 


i TT pyap-—2 1 fear (26) 


|2|/z a ; 


Writing |x|/r=&, one obtains (see Appendix 5): 


-[4 4 i? ( it Be) Vie — & pen (62 — 2k)? V1— &] + 


ag abo 7eS Be 


+[(1—4h) Bo+4h] F (resin V1 — 2, aresin Ve 2 wats \e 


Bo 


—4k (1 — ps) F # ( aresin y/1— (60° _ (4). arcsin 


i) 
5 at fo) £ 2 (aresin ae (*). aresin res F: | “7 (27) 


for By<&<k, and 


Ae 1 pt 4k? (1+ kh’) Bo +8k*|E (aresin Vi —p3 — Bs, aresin eee a i) a) - 
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2__ 9 R2)2//] — & 
0, (€)=— {-@ + 


~ 9g (Bo) EVE — B 
nas eae gers sap ee 
+ (1-44) Bo +4] F (aresin i= Beata hoon a 
; poeple 
—7ali—4 (1+ 2K) Bo+8 kt] B (sven ree aresin a) di 
bine Ok (Bo) (28) 


40° 9 (Bo) 
for k<&<1. There 
9s (Bo) = (1 — 4 #*) 2+ 4 k*] K (V1 — B3) — (1/B5) [Bo — 4 &* (1 + b*) 86+ 8 k*] E (V1 — Bo) 


and F and £ are elliptic integrals of the first and the second kind, respectively. 
Further 
oz(§)=—q for 0<&<f, (29) 


and oz (&)=0 foray Lene. (30) 


For By<1, §<1 one obtains 


é |x| 
(0 ee?) Le. Gz @)=0( jaa) 


in agreement with the well-known static solution, if 2v9t is identified with the 
width of the uniformly loaded strip. 

The stress o, at the surface z=0 has thus been completely determined by 
the expressions (27)-(30). An alternative expression, suitable for certain uses, 


is given by: 
sp-3) -5V =) : “a I 
_4O 2. 2k? at x xv x Ok? 
| «| / 7 r\3 
T(z) (5-1) (.- $5) 
' 


where ['(t)=T,+y+T,(t) is the path of integration shown in Fig. 5. The 
proper branches of the integrand are found by making indentations below the 
points t=|z| and t=|2|/k. 

It can be shown if y is sufficiently small that the expression (31) expresses 
oz for all values of x, except for the neighbourhood of |x|=£,7 (Appendix 6). 


2 


dt, (31) 


Problem D 


A pressure qd, is acting on the infinite strip |a|<v)t of the surface z=0 of a 
semi-infinite solid z>0. The displacement, w,, at the surface z=0 is zero for 
|z|>v,t. Find the displacement u, at the surface z=0 for |a|<vyt. 
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Lom t> % idmT (5 a = 
(o) (2) 
a esis 
Ret g Ret 
tT, @) TY) 


Fig. 5. Path of integration in the t-plane. 


The displacement u, at the surface z=0 


We may use the expression (31) for the stress o, on the surface. Thus the 
problem is no longer a mixed boundary value problem and the solution will be 
completely analogous to the solution of Problem B. We need to calculate the 
Laplace transform of o,. 

Supposing y (see Fig. 5) to be a semi-circle with radius e-|x|, we get 


ee hel 


4 3 4 
ae | fA (yar rer, 


where { (yee : <3 
Tc 2 
faye 
h(y=1-0(r-| lee j2|) +0 (7-2 é 1) 
0 Bo 
and T= f(x)dr. 
Y 


Writing aed e ?*g,dt, we obtain for real, positive values of p: 


« Tt = 
dys. Ren | € “fren ‘\dx' e)avde—Rep | « “Pt Idt= 
= — Re [ ee Aa h(x)dr—Rep | « PF Idt= 
i) 6 . 
= —Re Takienegl a?" ¥ (z ar—Re p | ¢ PP Let, 
T' (co) y 0 


where (co) =I',+y+I, (ce) is the path of integration shown in Fig. 6. 
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iImT 


Fig, 6. Path of integration in the t-plane. 


Thus, for «<1: 


Ay 


HG s.- — Re | et (tr) dt +e ?!*l/> Re [ far- 
T'(00) Y 
—Re p | Hd pe “td +0 (p2?et)= 
|2|/Bo 


= —Re | e ’*#(t)dt +O (pa et), 
T' (co) 


where O indicates “‘the order of’. As the inverse Laplace transform of the last 
term is 


6 (t) O (x? e*) 
this term corresponds to a load 
=0 (4k Q|2x|-2*)-6(z) 


which certainly has a very small influence on the displacement u, on the surface 
for |«|<f,t, and this term can thus be neglected in the following. Moreover, 
in the limit e—+0, one obtains the exact Laplace transform of the stress (c,)2-0 
corresponding to the conditions in Problem D: 


2 4 kt 
s.=p | ©" o,de=— are [ e fizkds. 
4 s 


T' (0) 


In analogy with the solution of Problem B, we wish to invert Fourier 


integrals of the kind shown in Eq. (11). Therefore we first deform the path of 
integration: 


4 4 
S,= — een Tire. ers 
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Then, after putting t=«|2x|/(ip), we obtain 


p 2 p 
ag ps —2|/(p'+24)(2 +2) 


Vp? + a? Vp? + of BR) 


S,=—- 


cos a|a|da. (32) 


The analogy with the solution of Problem B is then quite obvious. It is 


merely to replace the expression (5) by the expression (32), i.e. to replace the 
quantity 


ee oe 
2ak’ och p*+ a* 


occurring in Eqs. (13) and (14) by the quantity 


2 2 
p - ay aoe 
PQp (Ate) a” V (p® +00) (B+ at) 
0c Vp? + 02 V(p? + a? 
We obtain for z=0: 


—ioo 


p w= 


pitieg 


0c | re iy 


After deforming the path of integration to the path shown in Fig. 6, we get 


Q enya a | ae aoe 
pram Srey? | U(-G,) geet a 
x 
ral 
Ulr-—-— 
eu, ( 
and thus ial @ ice EOE : 


After two integrations we obtain 


sa’ i Ss a ( =I 
p= Ue: (35) 
0c Bo! Be Bo 
The shape of the surface displacement is thus a semi-ellips 
2 2,424 42 
bcc QO Ca oO. as u,>0. (36) 
Bot” Q pot 
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If By><1, and ~yt=a, one obtains, after making use of the relations 


Q (Bo) < %o 
pe 2k) 


when fy 90 


R= Lae . where y = Poisson’s ratio 
2(1—») 
0 c= itd : where E=Young’s modulus 
(l+yv) (1-27) 
2 u2 
the following equation a 1 2 1 


2 
J2a—r) Bal 


in agreement with the well-known static solution (cf. for instance Sneddon, 
1946). 


The solution of Problem D provides also a verification that the solution (23), 
(24) satisfies Problem C. 


The stresses in the semi-infinite solid 
Because of the conditions of plane strain, one obtains 
oy= (1-2) ockd (37) 
Ca¥ Ode 


5 (Re cad. (38) 


The state of stress is thus completely determined by the following three 
quantities: 


Ode Og tee 


Let us start to determine ¢. We obtain in analogy with Eq. (11): 


2 
oo ( P 2) ,-2Vp* +e? +il2 
+00 : (Fatale ZV p? +a? +i\rla 
O= f Ae-VP+2 cog REP Sp pa EE 2 [ 2h 

0 Ca 


~c0 v Vp +o V(p? + a? 63) 


Putting «=px, one obtains 


1 ae gt 
2h? Q : (sat) sa eae: 
Decgitis 2 Re [ dx. 
QCa : Vi+x? V+ x? p28 


The path of integration can be deformed to the straight line from x=0 to 
176 


ARKIV FOR FyYSIK. Bd 18 nr 10 


_”, 


arctg m 
ext 


Fig. 7. Path of integration in the T-plane. 


and that branch of the hyperbola 
(Re x“)? (Im x)? 1 


2" x Re 


for which Rex>0 and Imx>0. (See Fig. 7). 
The real part of the integral along the straight line is zero and thus the 
path of integration is reduced to the hyperbola-branch. If then the substitution 


6=zV1+x7—i|2|x 


is introduced, the path of integration in the 6-plane becomes the positive real 
axis from 0=V2?+z2=R to 0=00. 


One obtains ee A) 
V6? — R? 
and thus 
oo 1 2 
As _2#Q egg | ae” | 
‘Re | U1 eae Ly (Lat 3d Bey sme atte 
0 
Tt iia 
2k°Q dé za Lae af 
then $B [alia Lg. 0 
a. ary: V68— B® LV (1 + 2 BO) ee 
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> 


As 


ET 


Z 


Fig. 8. Coordinate ¢ and % in the vicinity of the point |z|= By 7, z= 0. 


Now we investigate the variation of the magnitude of ¢ at a small distance 
et from the crack-tip with the angular location as given by the angle #, where 


|z|=Byt+et cos #, =et sin J. (40) 


See Fig. 8. 

Insertion of (40) in Eq. (39) gives, after higher order terms of ¢ have been 
neglected, and the substitution w=6/t has been introduced, the following ex- 
pression for e<]1; 

Tn me 

2k? du BR 2k 

a Re Ve B ca eae ge ee Ee 73: 
Bote cos o ‘i Bo iow +25 ut cos P+ 215 wu uVe=fasin o| 


It is evident that the main contribution to this integral is given by values of 
wu in the neighbourhood of w=1, and thus we can write 


pr 


1 1 
qe 289i 35 | du sa: 
eca V1—£fR Bo |1—u (1-25 cos 9-24 5 VI= Bi sin 8) 
ri Bo Bo 


After integration we obtain 


Q 2k 2U—f » 1 


gp © e 
eee BoV1— Bo le cos P+ 2% 


5/1 sn | 


Q 2h- ea be |/ Vi=Bi sin? 0 + cos 9 Bo sin? 9+ cos 3 
d th $ 
+e ne oD oe ca Bo V1— Vl—fR 1— £§ sin? 3 (41) 
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It is to be observed that when 6,0 Eq. (41) changes to 


2 Bo oo a 
or 2 0c ( ens 22 pe gst 
and thus, because of Eq. (38) 


mee to 2° cos eet (42) 


in agreement with the static solution (Sneddon 1946). 
o, and tz, are determined in a similar way as that for ¢. One obtains 


2khQ 2h ke 


Pp 


S.= - 


cos axda 


+o 2 2 eae wii hammer a Pig 
(Fare e-2V p+? ~ of Vp? +o? p Be te #VP +e 
Vp? + a? V(p? + a? 63) 


+0 2 

Bee = (eased Be Sai 

Bee ee 2b epic BEA [e-2Ve+e — ¢22V% +a] sin arda, 
p V(p? + a Bo) 


—-0 


where S, and T,, are the Laplace transforms of og, and T;,, respectively. S, and 
T,. contain two parts, one irrotational (characterized by e-#VP*+a*) and one equi- 
voluminal (characterized by e~*)?"/"*+@"), After the substitution «=p, the path 
of integration is deformed to a straight line along the imaginary axis and a 
hyperbola-branch in analogy with the case when ¢ was determined. After the 
substitution 


= 


“i + 4? —i|2|x 


in the equivoluminal part, the path of integration in the §-plane becomes as 
before an infinite segment of the positive real axis. One obtains 


t 1 ae 
—— 1 
dé (53 
4H 9 Re | =| + 


Vo? — BL V(1+ x? fo)*J, VO eto 


R 


(43) 


. A( +2) 
dé 2h 7 Ri 
seat a = Re ae Bids 
Ye-=°" J 


R/k 
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and 


4 1 
“Vite (s+) 
2h? 
tn -4k*Q- Re | - = 3 . | 2 Ra + 
|| V0°= RP iV (1+ x? 63)? V 6? —R22+i6|2| 


R 


t ie % al 1 _) ) 
x"? (|— +x 
—-4k*Q Re i mit 
‘Te 4 | Ye R? iVa+2 1+? pa) Vor-2 2+10 | 


x= a 


(44) 


It may be observed that no head-wave appears. This fact has a simple phys- 
ical explanation: The surface z=0 is here prevented to move in the z-direction 
when |z|>f,7, and thus also rotation of the elements at this part of the sur- 


face is prevented. 


In the vicinity of the point |z|=/,7,z=0, Eqs. (48) and (44) reduce to 


,& Bo (2? — p3) |/ V1 fb sin* 0 — Bo sin? b+ cos 3 
“Bi 2V1—p Fe 1— f§ sin? 3 


+218 Ve R vi = (Bo/k') ane cos 8 
1 — (B6/k*) sin® 3 


. x RQ Bo 2? — pt AR sin? fe ponds 
eit [2] Banh ces 1 — f§ sin? 0 


_ = 65/2) sin? — tae sin? } — cos # 
— (Bo/k*) sin? 9 d 


where ¢ and # are given by Eqs. (40). 
When £,—0, Eqs. (45) and (46) reduce to 


Bb 8 ee 28 
0. to / cos 5 1+ sin* ?— sin 3): e<l 


and Tr2> = 05/8 sin i 5 (cos? 5 —— sin? 0), ex) 


in agreement with the static solution (Sneddon 1946). 
Eqs. (39), (43) and (44) determine the state of stress completely. 


180 


(45) 


(46) 


(47) 
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Discussion of the propagation of the brittle crack 


The velocity of a brittle crack 


By superposing a state of plane strain characterized by o,= qo, 0, = 0, = (1— 2k’) q 
to the semi-infinite solid considered in the preceeding sections, ah then putting 
two such semi-infinite solids in contact with each other, the solution (35) in- 
ae the propagation of an elliptic crack. However, the propagation velocity, 

=f ca, of the crack-tip as given by Eq. (35) presupposes that the two semi- 
infinite solids are joined together by a joint resistant to any finite stress if the 
velocity v9<c,, where c,=k,¢g is the velocity of surface (Rayleigh) waves. This 
follows from Eq. (27), which expresses the stress in the joint (after g, has been 
added to the right member) because the term 


1 = fees a2 
222 SSS 3 2 2 2 2)2 2 
iVe_ pith (1 — Bo) Ve? — & — (68-2)? V1— 2] (49) 
becomes infinite for =f), i.e. at the crack-tip, and further it is positive in the 
neighbourhood of =f, if By< ks. 

On the other hand, the term (49) is negative in the neighbourhood of =, 
if k,<f,<k. This indicates that, if the propagation velocity of the crack, v. 
has some value between c, and c,, we must presuppose that the joint cannot 
resist any tensile stress at all. 

Thus, if the joint has some finite strength, the velocity of the crack-tip will be 
precisely v.=c,;. This could perhaps have been predicted intuitively. However, 
we have arrived at this conclusion by tacitly assuming that the material of the 
two semi-infinite solids possessed an infinite strength. If we now suppose that 
this material and the joint are of equal (finite) strength, we are obviously dealing 
with the original problem, the propagation of a two-dimensional brittle crack in 
an infinite elastic solid. 

In this case the correspondence between Problem A (Problems C and D) and 
the original problem is complete only if the stress state in the neighbourhood 
of |x|=f,7 is such that fracturing occurs only in the plane z=0. In the ab- 
sence of a reliable fracture criterion for a brittle material, we choose to in- 
vestigate the behaviour of the dilatation ¢ and the normal stress o for a hypo- 
thetic fracture direction, defined by the angle #, given by Eqs. (40). og» is given 
by the relation 


(Sleep ee aes cos s29— lly, sin 20. 


From Eqs. (38), (41), (45) and (46) ¢ and oy» can be determined in the vi- 
cinity of the crack-tip. One obtains for e<1 


2k Bo Bo. 0 50 
$e peg Stat Vb F, (Bo, 8) (50) 
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(1-2) oe ) 
von 2V1—R zi or Fy (Bo ®) 


4k VE BBV 1 — BS Fy (Bo/k, °) _2—p)| - (51) 


ai @P—R) Fy (By 8) 


lace 
Bo (1— k*) 


_ 2k YV1—Bo V1—fo aeae F, (Bo, 0) — Fe (Bo/k, a 

Bo (1 — k*) Fy (Bo; 8 : 
uae Fs (fy 0) | OB eee 
and F, (By 8 vad V1— Bo sin? 8 — cos o 


1—£§ sin? 3 


By the aid of Eqs. (50) and (51) it can be shown that ¢ and oy both are 
maximum for #=0 if B,<k;. If a certain magnitude of ¢ or os were decisive 
for the occurrence of fracture, then the fracture would follow the plane z=0. 
However, a real material shows anisotropy as regards possible crack-surface di- 
rections (the crack-surface will follow lattice planes) and inhomogeneities as re- 
gards the local “atomic’’ strength. Thus, in a real material, the crack would 
no longer follow exactly the plane z=0, but it would oscillate across this plane. 
The tendency for such oscillations seems to be rather pronounced because the 
angular maximum of ¢ and gz is rather flat (6 and o, behave approximately 
like cos 9/2). A few calculations based on specific values of k and f, indicate 
that this maximum will be somewhat flattened when f, increases. 

If the crack oscillates across the plane z=0, the component of the crack-tip 
velocity in the |2|-direction will be somewhat lower than c,. 


The shape of the running crack 


Eq. (36) shows that the crack is elliptic. The ratio between the semi-minor 
and semi-major axes is 


Q (Bo) 
0 ¢4 Bo 


A calculation of this ratio in the case of a material with Poisson’s ratio =2 
has been carried out. The result is approximately 


q 
1.837 7 for - p)p= 


q 
1.04 Fr for “BJS ky. 


The running crack is thus more eccentric than a static crack subjected to 
the same stress. 


182 


ARKIV FOR FYSIK. Bd 18 nr 10 


idms 


Fig. 9. Path of integration in the s-plane. 


APPENDIX 1 
Calculation of wo (r) 
From Eq. (20) we obtain 
1 for Bo<r<k* 


QO (r)= gg ee ey es we 
(r) AN tk Vl —r) (r—k) +2 (r—-2 BP ise Tae BE 
Ami 4 V(1—r) (7-H) 4 (r— 2)? 


1 1 
© (s 1 ds 1 
eye Ae bag wade Seema 
Bs Bs 


a — 


opts cia Vi—s) GB) 
r={ SS 6 2k*)?+i-44 V(1 —s) (s—#) 


Wg, 
8—r  (g—2k)? 1-4 V(1—s8) (s—k*) 


where 


The integrand has the branch-points s=1, s=k*, s=0, s=ki, s=ry+in, 
8=r,—ir,, where the four last values of s are zeros to the equation 


(r —2 k?)*+ 16 k8 (1 —1) (r—k*) =0. 
It can be shown that k, is the ratio between the propagation velocities of sur- 
face (Rayleigh) waves and irrotational waves, and that ry<0. 


By choosing a path of integration shown in Fig. 9, one easily verifies that 
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(7-22)? 4° V(L—r) (k?—7) 


2 
27% In (ey for Bo<r<k 
I= ; 
2\4 6 ae aoe 
222 In V(r— 2") wel. us rhea) for k?<r<l. 
r (r-- ks) 
i tes PERS 
in ( hoa i +: fis ks) : for Be<r<k* 
r—Bo (r—2k*)?—4k3 V(1—1) (kh? —1) 
Then @ (r)= _ fm * 
ag eee ia kts : for ?<r<l. 
Fr Bo Vir—2k?)* + 16K (1 — 1) (k?— 1) 
APPENDIX 2 
3) 2 b (s) 
Calculation o 
“ie 1 [em V a? (s) + 8 (s) 
By 
idm s 


hs 
N 


‘Re s 


Fig. 10. Path of integration in the s-plane. 


This integral occurs in i (19). It can be split up into the two integrals 


D - eae (1 —s) (s — Bo) - (s —2 k*)? 


(s—r) [($—1)* + rz] 


ds, 


I aE s)(s— Bo)"(1—8) 9, 


(s—r ) [(s— 15)? oa 


Yo and r, are defined in ye a” 
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idms 


Res 


Fig. 11. Path of integration in the s-plane. 


By choosing a path of integration shown in Fig. 10 for J, and in Fig. 11 for 
I,, the two integrals are easily calculated. It is to be observed that the poles 
$=?) tr, can be disregarded, because the corresponding residues cancel when 
adding J, and J,. 

One obtains 
(M —r) for po<r<k* 

ti yet ee) ee) V(r — k’) (r — Bo) ‘eee 
(M-—r)t+a-4k (2 +168 dr) (r—B) or <r<l, 


1,+1,= 


where M is a constant. 


APPENDIX 3 
(r—2 k*)? 
Demonstration that g(r) =Q,° : sq 8 not a solution to Pro- 
(r — ks) V(1—r) (r— Bo) 


blem C. 
This demonstration is completely analogous to the calculation of uw, in Pro- 


blem D. One obtains 


x 


Z 


Fig. 12. Surface position at a fixed time, corresponding to the false solution when f, <,. 


185 


K. B. BROBERG, The propagation of a brittle crack 


Fig. 13. Surface position at a fixed time, corresponding to the false solution when f,>k,. 


|x| 


ee a Ul 
1 Q; 0(: i) 1 ( i) for By<k, 


Peck y/ | Bay eealal ee a ) 
eu, kK as an 
eo | 
: u (rll) 
3. Oe Bo 
oe aay aE) 
* ox? 


The solution for By<k, clearly does not satisfy Problem C. As for the solu- 
tion for B)>k;, after two integrations one obtains 


for By>ks. 


Uz 


SS ap All, 1 Aye 
_ 1 Bo NBs Qiks (--2)-1n Bo Bo Bo} _ 
Bo — ks ks yi-# 7—lel) 
Bo\ ke 


- — =... 
Bo V(B5— ke) 
k, || / bi y/ n° 
t+ -VY1l-j3/?-s 
- (+4) in | —B2 6 ai, @rkelel 4, _& 


1-4 (e+ ZI) “VB BP BB 
Bo ks 


In Figs. 12 and 13 the surface position at a fixed time is sketched. 
It is seen that the solution for By<k; does not satisfy the requirement u,=0 
for |x|>f,t (although @?u,/ér?=0 for |x|>f,7 with the exception |2|=k, 7) 
and that the solution for 8,>k, does not satisfy the requirement that u, should 

be finite (although u,=0 for |x|> 7). 
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The occurrence of the false solution to Problem C seems to be a result of the 
change of boundary condition from 


2 
=0 for |z|>f,r to oa 


a =0 for |x|> Bor. 


APPENDIX 4 
Calculation of Q 


The quantity Q occurring in Eqs. (23) and (24) is to be determined from the 
equation 


Q- 'Ve~ Be Out eae 


oe’ 14 


Thus Q= ; 2 4 Po 
Vs—By 9 (s) 
Ty UL ep easel 
- | Vexil. Omen 


where ae Q is given by Eqs. (23) and (24). One obtains 


? 


ke 
flac fs, Dae Vs — By (8—2%)?—4 8 V(1—s) (Bs) 
eV rr soem sV(1—s) (s— B3)° 


0 0 


ds + 


Vs — By (¢—2%*)? oF 
Vs #8 V(1—s) (s — 63) 


kK 


1 ka 
(Vf) 0-28, gp { WE PIVE= sg 
» Vs (1—s) (s— Bs) 2 8Vs(s—Bs) 


Lim (s—2 k*)? [5 3 ( Bo(s—2k) res 
I Yip ae 1—s) ( oi V's (1—s) (s— 3) (1—s) (s— Bo) 


7 Vie—s_ r pV Bae 
—4}' ds+4k ds 
of iP (s— BP ae sVs(s—B3) — Bo) 


The second and the fourth integral in the last member are calculated in Appen- 
dix 5. The first integral, Z,, and the third integral, /;, are easily calculated by 
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~Ims 


iY 
N 


Fig. 14. Path of integration in the s-plane. 


9) Res 


contour integration, choosing a path of integration for J, as shown in Fig. 14, 


and for JZ, as shown in Fig. 15. The radius of the semi-circles at s=£ 


(r — >). When r—>f5 one obtains 


ae (s —2 k?)? Diam tee es 2 (Bo — 2 k*)? : 
fr eam BR RV pr Be) 
i 2V ik pi 
n- | eae ee at ae 
Thus we get 
Vs B 2 2, 4 
Ve—Bo 9 (8) 9, 9 Po- 4k Bot 4k! Ap —a5 
p m9 ds=n—2 RB) K (V1— p2) + 
Bo—4k (1+k*) p6+8k4 4 
2. — pe 
y BS (1 — A) germina 
SE x (Vi - (2) ) 38 16H 5 2 (/1- (6) 
2 2 
rit _ .Bo(1— Bo) 
an en Q= 4% 9(Bo) 
where g (Bo) = [(1 — 4.) 5 +4 4] K (V1 — f2) — 
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idms 


Re s 


Fig. 15. Path of integration in the s-plane. 


ea ere 
-ues-anlf @))e-an(li- 


APPENDIX 5 
Calculation of oa, (&) 


From Eqs. (23), (24) and (26) we get 


(r—2 ei)" dr—4 | == dr 
r Vr (1—r) (r— Ba)° a Vr (r— Boy 


where &°= 2?/z®. If &>k the last term should be dropped. Inserting p*=r, 
we obtain after some rearranging: 


2089 f Ld a tzse +a 20 | pte 

Qa = om 

_sit ap a ict VERB 
B) RAE — Bo) sje | RF Bo 


These integrals can all be reduced to elementary integrals and elliptic inte- 
grals of the first and the second kind. The result is given by Eqs. (27) and (28). 
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APPENDIX 6 
Alternative expression for az (&) 


By the aid of Eqs. (23), (24) and (26) it is easily shown that Eq. (31) expresses 
o, for all values of |z|>fot if the radius of the semi-circle y is sufficiently 
small. From the first equation in Appendix 4 we obtain 


1 
1 Vr—Bo a Q 
sera) Vr (r) dr 465 9 f 
Bo 
r Vr—B Ve 
1 r— Be aQ ( E 
Way _ = Atm os 
ait) Synch aise 2a 
Bote 


1 1 
1 ( v(r) 7 Qi thal Ve 
3 "7 dr %~ 33, * OR [ emar+o(%o) 


when e<1. The integral in the right member has been calculated in Appendix 4. 
We obtain 


2 


1 f(r) Q[(p3—2 k*)? 4 V1 — fi) ( — f3)] (i 
dr= q+ irre =). 
} = #8 V1 — fi Ve py °) 


Now we study the integral 
jar pReg | CREAN gs 
r Vr (L—7) (r—B3)° 


T(r) 


? 


where ['(r) is the path of integration shown in Fig. 16. 
Then obviously 


(Bo 2 k*)* = 4&8 V1 — Bo) (k= Ba), 
B3V1—B3Ve 


“3 | Parole) 


Bote 


I= -4ReQ | e'Sdg— 
0 


_Ql(gs- 2k 4 VO— A) (Bl _1 f ptr) Ve Ve 
;3 1-fiVe ais Taro (40)--a+0( ) 
Bite 
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tImr 


Rer 


Fig. 16. Path of integration in the r-plane. 


As the integrand of the integral defining J is purely imaginary for real positive 
values of r<f>, I must be independent of ¢. Thus I= —q. 
After inserting r=2?/t*, one finds that 


( 1 -5) -5 V(5- i) ton 
4k*Q Dk iee ger \a* a ik? 
— Re dt, 
x / 2 ey 3 
(5-1) (1-5) 


Dy 


Tie 


where [', and y are shown in Fig. 5. 
But [= —gq. Further 


esl EES 
D2 nw2) ned oe oy 
Re 2k 2 x x a2 k d 


T(z) 


where I’, (1) is the path of integration shown in Fig. 5. This follows from the 
fact that the integrand is purely imaginary for t>|z|/f,. Thus, if r>|2|/Bp, 
the right member of Eq. (31) equals I= —q , and Eq. (31) thus expresses g, 
for all values of |z|+ 7 if y is sufficiently small. 


Conclusions 


The investigations have shown that a two-dimensional crack will propagate in 
one plane in a linear elastic medium which is homogeneous and isotropic as 
regards stress-strain relations and fracturing characteristics. This conclusion is 
partly based on simple specific assumptions regarding the fracture criterion for 
such a material. The shape of the crack is elliptic, just as in the static case, 
although the eccentricity is somewhat larger than in the static case. The velocity 
of the crack-tip equals the velocity of surface (Rayleigh) waves. 
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In a real material, showing inhomogeneity and anisotropy as regards fracturing 
characteristics, the crack will oscillate during the propagation. Thus it may be 
expected that in a real material (even if it is perfectly brittle) the crack velocity 
will be somewhat lower than the velocity of surface waves. 


(Royal Institute of Technology, Stockholm.) 
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